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Abstract. We propose a compact and explicit expression for the solutions of 
the complex Toda chains related to the classical series of simple Lie algebras 
0. The solutions are parametrized by a minimal set of scattering data for the 
corresponding Lax matrix. They are expressed as sums over the weight systems 
of the fundamental representations of fj and are explicitly covariant under the 
corresponding Weyl group action. In deriving these results we start from the 
Moser formula for the A r series and obtain the results for the other classical 
series of Lie algebras by imposing appropriate involutions on the scattering data. 
Thus we also show how Moser's solution goes into the one of Olshanetsky and 
Perelomov. The results for the large-time asymptotics of the A r — CTC solutions 
are extended to the other classical series B r -D r . We exhibit also some 'irregular' 
solutions for the T>2n+i algebras whose asymptotic regimes at t — > ±oo are 
qualitatively different. Interesting examples of bounded and periodic solutions 
are presented and the relations between the solutions for the algebras D4, B3 
and G2 are analyzed. 
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1. Introduction 

The famous Toda chain model § §, |, f§ |, | @, |, §, |l[ §J) was initially 
introduced in order to study the nearest neighbor interactions in atomic chains. Soon 
it was shown that it also possesses interesting mathematical properties and that to 
each simple Lie algebra q one can relate a natural generalization to the Toda chain 
@ 1 1 1 0, |, || 0,1101, namely: 

§=£a k e-^, (1.1) 
fe=i 

where q — (qi, . . . ,q r ) is a vector in the root space E r of the algebra g of rank r and 
a>k, k — 1,2, ... ,r are the simple roots of g. One may view qk(t) as the coordinate of 
the fc-th particle and study the effect of their interaction. A number of results in this 



direction are known showing how the (real) Toda chain (RTC) (1.1) can be viewed 
as completely integrable Hamiltonian system and how it can be explicitly solved, see 

§JJ J J, @ J, 0, [1 & 



Recently it became known that generalizing the RTC model with g ~ sl(N) to 
'complex' particles (i.e., now qk(t) become complex- valued functions) allows one to 
describe the interactions in the ./V-soliton trains of the nonlinear Schrodinger equation 
in the adiabatic approximation. In this case each soliton behaves as separate entity 
('particle'); Reqk(t) describes its center of mass position and Img^i) determines its 
phase, for more details see 0, [If, g gc}. Thcsc facts draw 

our interest towards the 

study of the complex Toda chain (CTC) models when the dynamical variables cj'fc(t) 
become complex, while the time variable t stays real. 

It is well known that a large number of results for the RTC are trivially generalized 
to the CTC case by just making the corresponding parameters complex. These include 
the Lax pairs and the explicit solutions. However, since each 'complex' particle has 
two degrees of freedom their interaction becomes much more complex and qualitatively 
different as compared to the real case. In particular the set of asymptotic regimes for 
the CTC is much richer than the ones of RTC. In addition to the asymptotically 
free particle regime (the only one possible for RTC), CTC allows also for bound 
state regimes, mixed regimes, degenerate regimes, etc. These facts were reported in 

The present paper is a natural extension of fed]; it also contains the proofs and 
generalizations of the results in J2(J . 

There are several methods for solving the RTC which readily generalize also for 
CTC. The method in §, |l| [l| allows one to write down the solution as: 

(ffl,u k ) - (q(0)^k) = lnMe-^'K), (1-2) 

where tuk is the fc-th fundamental weight of q, is the highest weight vector in the 
fc-th fundamental representation R(uJk) and L(fi) is the Lax matrix evaluated at t = 0, 



see formula (2.1c) below. The right hand side of (1.2) has the obvious advantage of 
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being written in compact and invariant form. However it is difficult to extract from 
it explicitly parametrized solutions. 

Another well known approach to solving the Toda chain models was developed 
in JnJ. ft allows one to express the solution in terms of 2r constants. Starting 
from a comparatively simple expression for X\ = exp(— qi(t)) one then calculates 
Xk = e~xp(—qk(t)) as a determinant of a fc x fc matrix whose elements are determined 
by the derivatives of X\ , see [Q, [nj . One may also use a recurrent procedure to evaluate 
Xk. However this leads to rather complicated and difficult to analyze expressions. 

Our first aim in the present paper will be to analyze (1.2) and write it down in 
the form: 



W),w k ) =lnB fl;fe (t), 

B sM t ) = E ex P 

7er g (aj fc ) 



-2( 7 , C)t+ (£ 0l7 )l W* fc > (£ 7) 



(1.3a) 
(1.36) 



Here ( 7 , £) is the scalar product between the vector £ = (£i, £2, ■ ■ ■ , ( r ) and the 
weight 7 £ T g (u>k); Cs are eigenvalues of L(0) and we suppose that they satisfy 

Cfe 7^ Cj f° r k 7^ j. The components of the vectors £ and ipo = (tpoi, ^02, ■ ■ ■ , Vor) 
provide the 2r (complex) parameters directly related to the minimal set of scattering 
data T g of L(0), see formula ( p. 15 ) below; T g (ujk) is the set of weights of the k-th 
fundamental representation of q; W^(C,t) are i-independent functions which arc 
defined in Section 3 below. Thus the right hand side of (1.3a), (1.36) like the right 
hand side of (1.2) is invariant and at the same time is explicitly parametrized. This 
fact allows us to calculate explicitly the large time asymptotics of q(t): 

v ± t) = 0±+0 ± . 



lim (q(t) 

t — >-±oo 



(1.4) 



It is a well known fact |l5|, [l6| that v + — w Q (v~) — where Wq is the Weyl group 
element which maps the highest weight Wk of the fc-th fundamental representation 
R[oJk) °f m t° the corresponding lowest weight vector ui^ . We provide explicit 
expressions for ^ as functions of £ and also show that 

/3+(C) = /r K(C)), = M0o) - 0o. (1.5) 

In fact the solutions to the CTC related to a certain simple Lie algebra g may be 
derived in two ways. The first one is to cast the solution (1.2) in the form 



E 

7er s (wfc) 



((« fc |V|7)) 2 e- 2 ( f 'T>* 



(1.6) 



and then try to evaluate the matrix elements (wfelV^) in terms of the scattering 
data T g . This requires the explicit construction of V for each of the fundamental 
representations R(ujk) of g. 

The second possibility which will be used below, is to start with the well known 
solution of Moser for sl(N) with conveniently chosen N and impose on the 
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scattering data the involution that restricts it to q. Obviously both solutions must 
coincide. The proof of this fact is also one of the results in the present paper. 

In the next section we introduce the notations and analyze the properties of 
the fundamental representations of the classical series of simple Lie algebras and 
derive some useful relations between the matrix elements of the typical and the other 
fundamental representations. In Section 3 we prove formula ( |l.3c|) , ( 1-31 ) for each 
of the classical series A r - D r . In Section 4 we extend the results for the large-time 
asymptotics of the A r -CTC solutions also to the other classical series B r -D r . We 
exhibit also some 'irregular' solutions for the T>2n+i algebras whose asymptotic regimes 
at t — > ±oo are qualitatively different. We provide also some interesting examples of 
bounded and periodic solutions and analyze the relations between the solutions for 
the algebras D4, B3 and G2. 



2. Preliminaries 



In what follows we shall use the so-called 'symmetric' Lax representation for the CTC 
model ([H]): 

r 

L{t) = ( b * H k + a k {E ak + E- a „)) , (2.1a) 

fc=l 
r 

M(t)=J2a k (E ak -E_ ak ), (2.1b) 
fe=i 

where a k = ie"^)/ 2 and b k = \dq k /dt. For g ~ sl[N) we have a k = I e (9 fc +i~« fc )/ 2 . 
It is well known that to each root a from the root system A 9 C E r one can relate the 
element H a of the Cartan subalgebra t)Cg. Analogously, to q(t) — Re q(t)+i Im q(t) e 
1} there corresponds the vector q(t) = Re q(t) + i Im q(t) , whose real and imaginary 
parts are vectors in the root space E r . 

The integrals of motion in involution for the CTC model are provided by the 
eigenvalues, C, k — n k + in kl of L. The solutions of both the CTC and the RTC 
are determined by the scattering data for L(0). When the spectrum of L(0) is 
nondegenerate, i.e. Cfe 7^ Cj f° r & 7^ h then this scattering data consists of 

T ={Ci,...,CN,r 1: ...,r N }. (2.2) 

where are the first components of the corresponding eigenvectors v^ k ' of L(0) in 
the typical representation R(u>i) of g, N = dimi?(ui). If we combine all eigenvectors 
as columns of the matrix V then r k = V\ k and 

L(0)V = VZ, Z = diag(Ci,...,Ov)- (2.3) 
It is known that the eigenvectors of symmetric matrices L(0) with nondegenerate 
spectrum can always be normalized, i.e. following |ll|, [l2| we require that 

JV 

( V «y fc ))^]T(V; fe ) 2 = l, k = l,...,N; (2.4) 

s=l 
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and beside s V T — V 1 . Equation (2.4) determines r k up to a sign. 
From it follows that 

G{t) = e- 2i W* = Ve- 2Zt V-\ (2.5) 



and we can rewrite ( |l.2[ ) in the form (1.3a), (1.3?) with 

B 3;k (t) = (Hl^l7)) 2 e- 2( ^' (2.6) 

7er B (^ fc ) 

Thus our aim will be realized if we obtain explicit expressions for the matrix 
elements (wfe|y|7) of V in the fc-th fundamental representation in terms of T which is 
determined by the spectral data of £(0) in the typical representation R(uii) of g. 

The eigenvalues of L(0), and especially their real parts n k , which can be calculated 
directly from the initial conditions as it will become evident below, uniquely determine 
the asymptotic behavior of the solutions |2(J . We will extensively use this fact for the 
description of the different types of asymptotic behavior. 

The minimal set of scattering data for g ~ sl(N) is obtained from (2.2) by 
imposing on T the restrictions 53fc=i Cfc = an d 

N 

/c=i 

which follows from V T = V^ 1 . Therefore one may consider as Tx r the set 
Ta t = {(i, ■ ■ ■ ,(n;<Poi, ■ ■ ■ ,<Pon} , N = r + l, 



^k = lnr 2 k -^Y,\nrl (2.8) 

s=l 

Although the number of elements in T^ r is 2N (instead of 2r = 2N — 2) it is obvious 
that only 2r of them are independent. 

For the other classical series of simple Lie algebras the elements of T ( |2.2| ) satisfy 
symmetry relations, namely pojj : 

Cfc =-Cs, k = N+l-k, (2.9) 
r k r- k =e-^w k , (2.10) 

for k = 1,...,N where N is the dimension of the typical representation R(u>i) and 
the value of qi(t) at t = is determined through the normalization condition (2.7). 
The coefficients w k are time independent and are expressed in terms of Ci , ■ • ■ , G- as 
follows, see Appendix A. 

B r -series: N = 2r + 1. Note that in this case ( r+ i = 0, 
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and in addition to (2.1C) 



' r+l 



-91 (0) 



W r +1 



W, 



r 1 
1 = 11 272 • 

8=1 Ss 



(2.12) 



Inserting ( |2 . 1 0| ) ( |2 . 1 2| ) into (2.7) we obtain a quadratic equation for exp(— <Zi(0)), so 
it can be expressed in terms of T g . 
C r -series: N — 2r and 

w * = -5^ A 27*317* II 77^7^ (2-13) 



D^-series: N = 2r and 



fc-i 

n 



4C 2 



4C fc 2 " 



n 



fe+i 



4C 2 



(2.14) 



In the last two cases again exp (— g'i(O)) is determined from (2.7). The derivation of 



the solution for the D r series requires some additional efforts. The main problem here 
is related with the t reatm e nt of the spinor representations. 

The proof of ( 2.11 )- (|2.14 ) is based on the study of the properties of the 
corresponding matrices V and is given in the Appendix A. Then one easily finds 
that the set of parameters 

% = {Cl,-",Cr;¥>01>--->¥>0r}, <PQk = M r kAfe)> ( 2 - 15 ) 

uniquely determines T ( [2.2] ), which in turn provides the full set of eigenvalues and 
eigenvectors of L(0). 

Next we will need a number of details from the representation theory of the simple 
Lie algebras. In what follows by R 3 (u>) we will denote the representation of g with 
highest weight u; r fl (w) stands for the set of weights of R g (uj). Often when the choice 
for g is clear from the context, we will omit the subscript and will write simply T(lv) 
and R{uj). We will also need to introduce ordering not only in the root system A B but 
also in the weight system R a (w). To this end we will use a vector K in the root space 
E r such that (71 — 72, K) 7^ and (a — (3, K) ^0 for any two weights 71 7^ 72 G T fl (u) 
and roots a 7^ (3 € A fl . Without restrictions we can choose K, together with the 
vector — k = — ReC to be in the fundamental Weyl chamber so that (u> — 7, K) > 
for any 7 £ r s (cj). 

Let us now denote by 7^, k = l,...,N the set of weights of the typical 
representation F 3 (ui) of g, namely: 



Ik = e fc 



r+l 
0=1 



N 



for 9 



7fe 



for 1 < k < r 

for r + 1 < k < 2r 



N = 2r, 



(2.16a) 



(2.166) 
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for g — C r , D r , 

!ek for 1 < k < r 

forfc = r+l N = 2r + 1, (2.16c) 

-ej. for r + 2 < k < 2r + 1 

for g ~ B r ; in all formulae above k = N + 1 — k. The corresponding weight vectors 
specify an orthonormal basis in R g (uji) and will be denoted by |7fc). 

An important and well known tool to construct the fundamental representations 
of g is to make use of the exterior tensor products of R g (cui ) . Indeed, the orthonormal 
basis in A fe i? fl (w 1 ) consists of the weight vectors 

hi) = hi u i2,-,ik) = hh A 7f 3 A . . . A 7i fc ), I = {k,i2, ■ ■ -,ik}, (2-17) 



with 1 < i\ < i2 < ■ . • < ik < N; the weight corresponding to (2.17) is obviously: 

li =7h +J l2 + -.. + 7i fc - (2-18) 

Here and in what follows we shall use the one-to-one correspondence between the set 
of indices / and the corresponding weight 7/ and weight vector I7/). 

For q ~ A r all fundamental representations are in fact exterior powers of the 
typical one R(u>i): 

R(u h ) - A fe i?K), (2.19) 

for k — 1,2, ...,r = N — 1 and 7/ = 7^' where the upper index k means that 
^( k ) g r(wfc). We remind also other well known fact, namely that in A k R(uji) we have 

MV| 7 ! fe) ) = (71 A 72 A ... A lk \V\ ltl A 7i2 A . . . A74) = V ( / 2 • • • * ) , (2.20) 
where V \ 3l J . 2 " ' Jk [is the minor of the group element V € © determined by 

1*1 «2 ••• Ik J 

the intersection of the rows ji,j2, ■ ■■ ,jk with the columns i\,i2, ■ ■ ■ ,ik- Thus given a 
group element V € SL(r + 1) in the typical representation i?(wi) one can construct 
its image for each of the fundamental representations R(u>k), k = 1, 2, . . . , r. 

Let us now explain how this can be done for the other simple Lie algebras of 
the classical series. To this end we shall make use of the well known facts about 
the root systems [ pi| p5| of g and about the tensor products of their fundamental 
representations, see }26[ 

Let us now consider the B r series. Then we have: 

A fc = R(u k ), for k = 1, 2, . . . ,r - 1, 

A r i?(wi) =i?(2w r ), (2.21) 

where 

w/s = ei + e 2 + ... + e fe , w r = |(ei + e 2 + . . . + e r ) (2.22) 
Here w r is the highest weight of the spinor representation of B r . Therefore the relations 



(2.20 ) hold also for g ~ B r with fc = 1, 2, . . . , r - 1. 
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Other well known fact is that the symmetric tensor product of R(u> r ) is generically 
reducible and 

[r/4] 

R{oj r ) ® R(uj r ) = R{2uj r ) © E R{v r -M+i) © R(u> r -u), (2-23) 

5 i=l 

where luq = and R(luq) is the trivial 1-dimensional representation of g. 

We have now two possibilities to introduce a basis in R(2u> r ). The first one is to 
use ( 2.20| ) and ( 2.21| ) like above. The second possibility is to argue that R(uj r ) eg) R(uj r ) 

is spanned by: 

| 7 « ® ^) and -L (| 7 M ® 7 00) + | 7 M 7 « } ) ? (2 . 24) 



where 7 | r ^ € T(o;,) are weights of the spinor representation of B r . Obviously they 
have the form: 

^^E^-^E^, (2.25) 
fe=i fc=i 

where <7fc = ±1. The corresponding sets I (and J) now must be special in the sense 
that: a) r+1 ^ /; b) if 6 / then ik = N+l—ik £ I- In other words I does not contain 

(r) (r)\ 

pairs of 'conjugated' indices. To the weight vector | 7 ) <g> 7} ) there corresponds the 

(r) 

weight 27) which can be obtained from 2to r by a Weyl group transformation. This 

means that | 7 j (^jj^) belongs to the module R(2oj r ) in the right hand side of ( 2.23 ). 
Thus we find: 

1 2 ... r 



= (u r ®u r \V\^ r) ® 7 i r) ) - (K|y| 7 i r) )) 2 • (2.26) 

Next we choose g ~ C r series. It is well known that the exterior products of 
R{lj±) generically are reducible, namely: 

A fe i?M = R(w k ) © A fe - 2 i?M, 

w* = ei + e 2 + ... + e fe , for fe = 2, . . . , r, (2.27) 

but the re latio ns fl2.17D , ( |2.18| ) and Q2.20[ ) hold also for g ~ C r with fc = 2, . . . , r. 
Equation (2.27) reflects the existence of a nontrivial invariant subspace in A 2 i?(<x>i), 
see @: 

JV iV 

|c) = E(- 1 ) fe+1 |^ A ^)= E ^17* A 7m), (2.28) 

fc— 1 fc,m— 1 



where S is the matrix entering into the definition of the symplectic group (see (A..5) 
below); namely 

XeSp(2r) «-> SX T S- 1 =X~ 1 . (2.29) 
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It is easy to check that due to ( 2.2S ) we have X\c) = \c) for any element X S Sp(2r). 
Thus we establish that c determines the one-dimensional invariant subspace in 
A 2 i?M = (C|c))ffiiZ(w2). 

Let us now analyze the weights in the weight systems r(cjfe) and their 
multiplicities. The highest weight e± + . . .+efc has multiplicity 1; the corresponding set 
of indices is I — {1, 2, . . . , fc}. Let us consider next the weight jn^ = e\ + . . . + Ck-2] 
to it there corresponds each of the following sets of indices I = {1, 2, . . . , fc — 2,p,p}, 
k — 1 < p < r. Therefore to 7^ there corresponds the subspace V(7m) C A fe i?(wi) 
which is spanned by the vectors I71 A . . . A jk-2 A j p A 7p) and has dimension 
dimV / (7( 1 j) = r — k + 2. At the same time the Freudenthal formula shows that 
the multiplicity of 7m = e\ + . . . + &k-2 in R(u>k) is r — fc + 1. This difference 
is due to the fact that in V(7(i)) there exist a one-dimensional invariant subspace 
determined by I71 A . . . A "fk-2 A c). The same argument can be applied to each 
of the weights 7^ = 1117(1) where w is an element of the Weyl group. Indeed, 
it is known that the Weyl group preserves: i) the lengths of the weights, i.e. 
(7(i) > 7(i)) = (7(i) '7(i)) = k — 2; ii) the multiplicities of the weights. In fact the Weyl 
group is isomorphic to the group of permutations 6>2r of the indices {1, 2, . . . , 2, 1}; 
therefore instead of looking at the transformed weight 7^ we may consider the 
corresponding set of indices which can be obtained from / by applying a specific 
element of 6>2 r - This analysis can be continued also by considering weights of length 
k — 4, e.g. 7(2) = &x + • ■ ■ e^_4. Skipping the details we formulate the result, namely 

A fc R{u x ) = R{uj k ) 8 (C|c» A (A fc - 2 i?(^i)). (2.30) 
It remains only to note that from ( 2.20[ ), J2.30 ) and from V\c) — \c) there follows 

(iv k \V\ lh A...A lik )=V 1 ' ■■• L 



ii ...ik 

H Wfe |y|7j fe) )+PM^|cA7/') 

= («>k\V\^% (2.31) 

where p is some constant, jj k ^ 6 R{wk) and 7/' 6 A fe_2 R(u>i) . In the last line we used 
also the fact that the representations R(u>k) and A k ~ 2 R(oji) span mutually orthogonal 
subspaces of A k R(uji). 

Finally let g ~ D r series. Then we have: 

A fc i?(^i) =R(u) k ), forfc = l,2,...,r-2 (2.32a) 

A r_1 R(u>i) = R{uj r -i + uj r ), (2.326) 

A r R(ui) = R{2uj r ) 8 i?(2w r _i), (2.32c) 



where 



^fe = ei + e 2 + . . . + efe, for fc = 1, 2, . . . , r - 2, (2.33a) 
w r -i = ^(ei + e 2 + . . . + e r _i — e r ), (2.336) 
w r = |(ei + e 2 + . . . + e r _i + e r ) (2.33c) 



The complex Toda chains and the simple Lie algebras 10 



From ( 2.32a ) and ( 2.33q ) we find that the relations ( 2.20 ) hold also for g ~ D r with 
k = 1, 2, . . . , r — 2. Let us now analyze the spinor representations R{u) r -\) and R{u r ) 
of D r . It is known that: 

[(r-l)/2] 

R(u r ) <8> R(u r -i) = R(u r -i + w r ) © R(u r -2i-i), (2.34a) 



»=i 



[r/4] 



R(cj r ) ® R(uj r ) = R(2u r ) © V i?( ). (2.346) 



i=l 



where <S> means the symmetrized tensor product and R(luo) is the trivial 1-dimensional 
representation of g. 

The basis in R(u> r ) ® i?(w r _i) is determined by |7j (g) Tz^" 1 ^) where 7j S r(o; r ) 
and 7}T -1 ^ 6 r(w r _i). The sets of indices I = . . . ,i r } and /' = , . . ,i' r } are 
related to the weights and "fp -1 " 1 by 



- | 2 ^ = iEt%- ^ = S E <™ = I E 7fc , (2-35) 
fe=i fe=i fe=i fe=i 



7/- - 2 
where cr^ and o~' k take values ±1 and 

r r 

JJ o-fe = 1 and Yl = L ( 2 - 36 ) 

fc=i fc=i 

Let us now consider a pair of weights 7}^ and 7}^ such that 07. = Oj. for all but 

one value of k, i.e. I C\ I' = {j±, . . . , j r -i}- Then the vectors 7j + 7}^ X ^ and 
w r + uv_i have the same length and one can check that they are related by a Weyl 
group transformation. Therefore {jj^ ® 7^ belongs to the module R(u r -i + uj r ) 



( ^.34o| ) and 



= {^V^-^AV^P). (2.37a) 

Finally, the spinor representation R(u> r ) is considered quite analogously to the one of 
the B r case. There is only one additional fact that should be taken care of: now the 
weights 2uv and 2w r _i have the same length but are related by an outer automorphism 
rather than by a Weyl group element. In order to separate the weights belonging to 
the modules R(2u) r ) and i?(2w r _i) we need projectors onto each of these invariant 
subspaces in A r R(u) r ) in ( ^.32c ). As we shall see below we will need to sort out only 
the weights of length r which are given by either 27} or by 27} . The projectors 
which will separate these two types of weights can be constructed by using (2.36) and 
the fact that Q k = —Cfc! the matrix elements of these projectors are given by: 

ft ^ fti,,.,i r = s( i± ot.'.'i) ■ (2-376) 
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Indeed, it is easy to check now that 

r+ (l — 1) r — (r) n r+ {r 

fill' = Ji'Ti =0, fJTj 
Thus we obtain the relation: 



7/ 



fi'Ti' 



= 7/< 



1/ 



1 2 

«1 «2 



(2 Wr |V|2 7 «) = (K|V| 7 j r) » 2 . 



(2.37c) 



(2.37d) 



From ( |2.37Q )-( 2.37a] ) we find the necessary expressions for (w r _i|y|7j 1 ') and 
(kV|V|7j ) through the minors of V, see ( A. 14 ). 



3. The solutions of the CTC revisited 

We will now analyze the structure of the solutions to the CTC for each of the classical 
series of Lie algebras. Our aim will be to write them down in the form (1.3a), ( 1.3^ 
and calculate the functions W^((^,j) for each of the classical series separately. 

3.1. The A r series. 

It is well known j|, [| ^2|, [33) that the solution s in this case can be expressed through 
the principle minors of the group element G(t) ( |2.5| ) . Using the Binet-Cauchy formula, 
or equivalently ( 2.20| ) and (AJ5) we obtain the functions: 

N 

A 1 (t) = J2rte- 2Ck \ (3-1) 

k=l 



Mt) = E r l r l ■ ■ ■ rle-^^-^W 2 ^, . . . , i k ), 

ii<i2<—<ik 

N 



A N (t) = ]Jr 2 s W 2 (l,2,...,N) 



-iVgi(O) 



(3.2) 



(3.3) 



which are proportional to the above-mentioned principle minors of G(t) . Her e 
r k are the scattering data of the Lax matrix L(0) introduced in (2.2)— (2.7) and 
Wi(C) = W(ii, i2i ■ ■ ■ ,ik) is the Vandermonde determinant: 



w I (o=w(i 1 ,i 2 ,...,i k )= n 2(c,-cp). 



(3.4) 



Note that each of the factors in the right hand side of (3.4) can be viewed as the scalar 
product (C)Ct) where a is an appropriately chosen root a € A fl . This and the fact 
that (£, Wo(a)) = (wo((),a) make it easy to introduce in a natural way the action of 
the Weyl group element u> on Wi(£), namely w : Wi(() — » Wi(wo(q). 
The solution of the corresponding CTC is then given by: 

(q(t)-<l(0)^k)=lnA k (t) (3.5) 
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g*(t) = gi(0) + ln 



A k (t) 



(3.6) 



with Aq = 1 and gi(0) denned by (3.3). 

For further convenience we introduce the functions 

B k (t)=e k ^A k (t), (3.7) 

and rewrite the solution to the CTC in the form: 

(0),u k )=ln.B k (t) (3.8) 

Introducing now th e se t of variables ( 2~8|) a nd ta king into account ( 2 . 1 6 o| ) we 
easily cast the solution ( |3.8| ) into the form ( 1.3tj ), ( |l.36| ) with 

B k (t)=B k (t) = Yl e X pfa), 7 P)w?(()(W(l,2,...,N))- 2k / N , (3.9) 



7? ) er(u, fe ) 



¥>(t) = ('flit), ■ ■ -,<PN{t)) ■ 



<Pk(t) = ~Kkt + (pQk, 



(3.10) 



where tpofc were introduced in 
with 



Another possibility is to use 0(t) — —2(t + cpo 



2 r z 
Vok = &ok-kiw k - — ]nW(l,...,N) = qi(0) +ln-^ 
N w k 



iv k 



k-i 

n 



i 



JV 

n 



i 



which gives: 

JV 

il<J2<...<ir 

with p k — e Vok . Then finally Bx r ;kif) can be rewritten as: 



7j"°er(u> h ) 



with 



wW((,y k )=w k , wW((,^)=W?(t)l[w s , k = 2,3 



sei 



Now it is more natural to consider 

TAr = {Cl> • • • ) CjV, ^01, ■ ■ ■ , VOJv}, 



N = r + 1 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



rather than (2.8), as the minimal set of scattering data for the A r algebra case. Indeed, 
the elements of fl3.15| ), like the ones of ( |2.8| ), naturally satisfy the identities: 

JV JV 

Ck = (C, e) = 0, Y, = (fl t) = 0, (3.16) 
fc=i fc=i 
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e = e\ + . . . + e^v, which is directly related to the fact, that all the roots of A r also 
satisfy (a,e) = 0. 



However (3.15) has also the advantage that imposing on its elements the symmetry 
condition (2.9) and 



k = N + l-k, 



(3.17) 



one can get the minimal set of scattering data for the B r and C r series. 

Indeed, let us c hoose g ~ A 2r and let us impose on 7A 2r (|2.9D, ( 3.17 ). First note, 
that from ( |2.9| ) and (3.11) with N = 2r + l we immediately find that Wk takes the form 
(2.11) for k = 1, . . . , r and ( 2.12| ) for fe = r + 1; besides one can check that wj, — w^. 



Next the condition fl3.17| ) with ip Q j, = <?i(0) + ln(r^/w fc ) leads immediately 



to the relation (2.10). Expressing from it qi(0) as Yo.{wk/ {r^k)) we easily get 
¥?o,fc = l n ( r fc/ r fc) f° r = 1j ■ ■ ■ i r an d t/5o,r+i = 0. Thus we showed that 7a 2t , with 
( p| ) and ( |3~17| ) reduces to (|1|) for the series B r . 

The same procedure applied to T^ 2r l provides ( 2.1 5| ) for the series C r . 



3.2. The B r 



The solution is obtained from the one for A 2r series by imposing the relations (2.9)- 
(2.12). Due to this only the first r of the functions Ak(t) in (3.5) are independent. From 
the analysis in Sect. 2 we see that only the expression for A r (t) requires additional 
special treatment||, see (2.26). Thus we find: 

A k (t)=e- k ^B k {t), for fc = l,...,r-l (3.18) 

Ar(t) =e- rqiW B 2 r (t). (3.19) 

where Bk(t) are of the form < \1.3b\ ), 



Bk(t) 



E 

7l fc) er( Wfc ) 



exp 



(3.20) 



with 



sei 



for k = 1, 



- 1 and 



<p(t) = -2{t + <p , (3.21) 



(3.22) 



s=l 



Here Wi(C,) is the Vandermonde determinant (3.4). We remind that to each weight 



7^ <G T(cjfc) there corresponds an ordered set of indices I = {ix< ... < ik} and that 



Wk are defined by ( 2.11 ), ( 2.12 ). Each set / uniquely defines jj . If the weight jf' has 



„( fc ) 



multiplicity greater than 1 then there exist several sets of indices related to 7} 

X This is related to the existence of the spinor representation. 



For 
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example, the weights e± + . . - + e k and e± + . . . + efe_i in T(LJ k ) have multiplicities 1; the 
corresponding sets of indices are {1, 2, .... fc} and {1, 2, . . . , k — 1, r + 1} respectively. 
The weight ei + . . . + e/c-2 however has multiplicity r — fc + 2; one may assign to it each 
of the following sets of indices {1, 2, . . . , A; — 2,p,p} where k — 1 < p < r. Analogously 
to the weight e\ + . . . + e k -3 and e\ + . . . + e^_4 with multiplicities r — k + 3 and 

^ j one can assign each of the sets {1, . . . , fc — 3,p, r+l,p} with fc — 2 < p < r 

and {1, . . . , fc — 4,pi,p2,pi,p2} with fc — 3 < p\ < P2 < r. In other words we find that 
the number of sets / corresponds to the number of weights provided each weight is 
counted as many times as its multiplicity. 

We note that B k (t) can also be written in the form: 

B k (t)= 2cosh(V(i), 7 f )iy«(C,7f ) )+ £ ^(CWf), ( 3 - 23 ) 

7^>0 /: 7 < fc) =0 

where the second sum runs over the sets of indices corresponding to the weight equal 
to zero. For the B r series these sets are of the form: 

{p 1 ,...,p s ,p 1 ,...,ps} for fc = 2s, 

and 

{pi, . . . 7 p Sl r + l,pi, . . . ,p s } for fc = 2s + 1. 
Formula ( |3.23 ) reflects the fact that T(uj k ) is symmetric in the sense that if 



7j G r(w/.) then — jj = Wq^^) G T(uj k ). One can also check that 

w^(li\ k) ) = w^(wo(0,Mii k) )) = v^ (fc) (C,7f ), 1 = {?*,... M- 

3.3. The C r series. 



Remark 1 When one imposes the symmetry condition (2_S), ( 2.1L ) on the A2 r -i 
CTC the corresponding system of equations is slightly different from jjl.j ). The 
difference consists only in the coefficient of the term e - ^'"''- 1 which comes out as 
a r /2 instead of a r . The extra 1/2 factor is easy to take into account and therefore for 
the C r -series the relation between qk(t) and B k (t) is slightly different, namely: 

q k {t)=\n^- + \\n2. (3.24) 
Ok-l(t) 2 



Now we insert ( [2.13 ) into (|Q|)-(3.2) to get the solutions for the C r series. Thus 
we get 

A k (t)=e- k ^B k (t), (3.25) 
for all fc = 1, 2, . . . , r where 

B k (t)= e X p(m,li 1 ..A k )WW({,'y il ... ik ), 0(t) = -2& + o , (3.26) 

ii<...<i k 

k 

W^(C,7k)=w k , W^(C, lll .., lk ) = W 2 (i 1 ,...,i k )l[w ls . (3.27) 

s=l 
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and 7i!...i fc is a weight in A k R(uii). From our analysis in Section 2 (see ( 2.3C ) each of 
the weigh ts Ji lt ...,i k — 7^ + ■ • • + 7i fc can be split into |7 H ,...,i fc ) = |7^) + p\c A 7/') 
(see ( 2.31 )) and we have 

1 ■"M = { U *|7| 7 f ) ). (3.28) 
^i ■ ■ ■ A fc j 

This means that the summation over all sets of indices i\ < i2 < ■ ■ ■ < ik in (|3.26|) 
reduces to the sum over the weights of T(ojk) only: 



{uk\V\~fi, ik ) = V 



Bk {t) = B k {t)= E ex P m),ri )w (k Hc,iT y ) 



(3.29) 



2,3, 



(3.30) 



s£l 



i.e. we cast the solution in the form (1.3a), (1.31). Besides, like for the B r -series again 



-7/ so we have: 



B k (t)= 2 cosh ((^(i), 



(fc) 
7i 



M/W(C,7f) 



E ^ (fe) (c,7l 



(3.31) 



7< fc) >0 



1-7} - 



We note that the sets of indices corresponding to the 7^ = are given by 
{pt, . . . ,PsiPi, ■ ■ ■ ,Ps\ f° r k — 2s and by the empty set for k = 2s + 1. 



3.4- The D r series. 

There are some differences in treating this case due to the fact that the Lax matrix L(0) 
( |A.8D is not a tri-diagonal one. Nevertheless Moser formula ( 3.1 )—( 3.S ) for N = 2r, 
together with the corresponding involution (|2.9|), (2.1C), (2.14) provides the solution 



to the D r . Due to the somewhat different structure of the eigenmatrix V we find that 



the first r — 1 functions A k (t) are given by (jOJ), (3.2) and only A r (t) must be replaced 

by & 



Ar(*)= E ^■■■^exp(-2(0 1 +... + CvW(/+..., l J 2 ^ 2 (H,...,z I -). (3.32) 

ii<...<i r 

Besides, the D r algebra s have two spinor representa tions w hich requires additional 
care. Note also that due to ( A. 14 ) the projector f£ ; ( 2.376|) enters in a natural way 
into A r . Thus in the right hand side only the terms related to the weights of R(2u> r ) 
give non- vanishing contributions. 

Let us introduce the variables fl2.15 ) and along with (3.7) for k = 1, . . . , r — 1 let 
us put B r (t) = e rqi ^ A r (t); then we can rewrite the solution for the D r -series in the 
form: 



(q(t)^k) = lnB fe (t), forfc= l,...,r-2, (3.33) 

w r _i +w r ) = lnfl r _i(t), (<f(i),2cj r ) = ln5 r (i). (3.34) 
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Our analysis in Section 2 shows that w k = and 

B k (t) =B fc (t), forfc = l,...,r-2, 

B r -l(t) = B r -l(t)B r (t), B r (t) = B r {tf , (3.35) 

with 

B k (t)= ^p(mA k) )w^(C,4 k) )> forfc = l,...,r, (3.36) 

7l fc) er( Wfe ) 

Here: 

W (1) (C,7fe)=^fe, ^) = -2Ct + ^o, (3.37) 

W W (C, 7^) = W?(£)l[w., k = l,...,r-2, (3.38) 

7 j r " a) ) = V^i • ■ • w r W!(0, a = 0,1, (3.39) 



and the set of indices / = {ii, . . . , i r } in (3.39) is such that it determines uniquely the 
weight 7^ r ^ G r(w r _ a ) in the corresponding spinor representation. In addition 

(gtt),w r _„)=lniB P _ (t). (3.40) 

The explicit solutions for RTC with the simplest choices for g with rank 2 in 
invariant form were proposed in the monograph |ll ] ; they coincide with the particular 
cases of the ones given above provided a proper identification of the variables is 
performed. Here we choose to provide as examples the solution to the D4 case and 
its relation to the solutions for B3 and G2. 

Example 1 Let g ~ D4 = so(8). Then it has three 8-dimensional representations: 
R(loi) and the two spinor ones R(u>3) and R(u>i). The representation R{u)2) is of 
dimension 28. We remind also that D4 has an outer automorphism v\ of order 3 
which interchanges the 8- dimensional representations; more precisely: 

V\ : OL\ — ► 0:3 — » Qf4 — * CL\\ V\Oi2 = Ct2i 

V\ : lui — ► lu 3 — > LU4 — ► u>i; V1CO2 — (3-41) 
The equations of motion for the Ti^-CTC have the form 

q 3tt = _e«»-» + e 94 ^ 3 + e- q3 - qi , q 4>tt = -e qi ~ qs + e- q3 ~ qi . (3.42) 
The solution is provided by 

B 2 (t) 



?i(t) = lnBi(t), ffi(t)=ln 



B 1 (ty 



B 3 (t)B 4 (t) B 4 (t) 
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where 

4 

B 1 {t)=2j2 w k cosh( Pk(t), (3.44a) 



fe=i 



B 2 (i) = 8 X! WjWj [(Cj - 0) 2 cosh(^ + y>j) + (0 + Ci) 2 cosh(ipi - <pj)] 

4 

-16^^ 2 C 2 , (3.446) 



l<3 



4 i r 

B 3 (t) = 2- 11 [] _-j W(l, 2, 3, 4) cosh 



Lp 1 + tp 2 + Lp 3 - *P4 



>y >2 v / o 

2, 4, 3) cosh ^ + " + ¥>4 + 3 , 4, 2) cosh ^ = y2 ± y3 + ^ 

+W(2, 3, 4, i) cosh ^^^i j , (3.44c) 

2 



4 i r 

B 4 (t) = 2- 11 J] _-j W(l, 2, 3, 4) cosh 



2, 4, 3) cosh ^ + ¥>4 + ^(1,3,4, 2) cosh ^ = ^ 2 ± y3 = ^ 4 

+W(2, 3, 4, 1) coBh ^-^-^ + ^ j , (3.44d) 



tp k (t) = -2( k t + ip ok . (3.44e) 

We remind also that 4 = 5, 3 = 6, 2 = 7, 1 = 8; W(i,j,k,l) are defined by (3.4) and 
in both the summation and the products, denoted above by i < j we mean that i and j 
take values from 1 to 4- 

It is well known that the automorphism V\ maps not only the fundamental weights 
ujk as in (3.41) but also the whole sets of weights, e.g. v\ : T(u>i) — ► T(uj 3 ) — > T(oj4). 
Using the definition of Bfe(i) and the properties of the scalar products (C, u i7^) = 
(ff 1 Ci t'*') we obtain 

Bi(t\C,<pb) = B3(t;vi 1 C,Vi 1 (p ) = ^(^iOi^o), 

B 2 (t; C, #>) - B2(t; «iC, vi<Po)- (3-45) 
These relations are compatible with Jl.3c| ), ( 1.36|) , i.e., 

(fl[t),W3, 4 ) = InB 3 ,4(t;C»W)) = (#(*)> vfV) = lnBifcvf^vf 1 ^). 

Example 2 Aferf we will derive the result for the B 3 starting from the D4 case. It is 
well known that B 3 can be obtained from D4 by imposing a symmetry condition with 
respect to the outer automorphism v 2 of D4 defined by: 

v 2 a k — ak, fork = 1,2; v 2 a 3 — a 4 , u 2 a4 = a 3 . (3.46) 
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The symmetry with respect to v 2 reflects on the scattering data of the D 4 -case by 

C 4 = 0, <A)4 = 0, or ^(t)=0. (3.47) 

To the end of this example, when referring to the variables related to D4 we will 
assume these conditions imposed and will denote this fact by additional "prime"; 
the corresponding variables for the ^3- case w e wil l denote by the same letters with 
additional "tilde". Then inserting ( 3.^ r i ) into ( 2. 14 ) we can write: 



W4. 



w' k = 2w k , k = 1, 2, 3; v 
In analogy with ( 3.4ty we have 

Ki^Cfo) = B 3 (t;v 2 (,v 2 <f ) 
and after imposing the v 2 -involution symmetry, namely, (g(t),ui k ) 



(3.48) 

(3.49) 

(q(t),v 2 uJk) for 

each k we have q' A = 0. Indeed using (3.4' ) one can find that B\(t; £, (po) — B' 3 (t; £, <^j) 
which leads to g 4 = 0. This way q[, q' 2 and q' 3 give us solutions of a system, equivalent 
for that 0/B3. More precisely after some rearrangements we find that 

B[(t) = 2B 1 (t), B' 2 (t) = 4B 2 (t), B 3 (t) = B' 4 (t) = 2 3 / 2 B 3 (t). (3.50) 

Comparing \3. 5l\) with ( 3.4$ ) we see that 

q' k (t)=q k (t)+]n2, fc = l,2,3; <&(t) = 0, (3.51) 

where 



$i(t)=lnBi(t), 



q 2 {t) = 



q 3 (t) -ln^ V ' 



(3.52) 



Bi{t)' " w S 2 (i) 

77ie condition ( 3.4ty , or equivalently, q' 4 {t) = when imposed onto the system 
of equations (3.4%) leads to a system for q' k {t), k = 1, 2, 3 slightly different from the 
B 3 -CTC. The difference is in the coefficient in front of a 3 e~^ q which comes out 
with an additional factor of 2. This factor is precisely cancelled if we go over to the 
variables <jk(t). 

Quite analogously but with more technicalities one can prove that the symmetry 
with respect to the outer automorphism v 2 of D r will reduce the D r -solution to the 
one for B,. 1; more precisely, using analogous notations as above, we can write: 

B' k (t) = 2 k B k (t), fork=l,...,r-2; B'^t) = B' r (t) = 2^/ 2 B r ^(t),(3.53) 

and 

«*(*) =&(*)+ In 2, fork = l,...,r-l, q' r {t) = 0, 



q k {t) =ln-~ 



B k {t) 



B k -x{ty 



for k = 1, . . . ,r- 2, 



B 2 r 

q r -i{t) = In -7T- 



(3.54) 



B r - 2 {t) 

In deriving ( 3.5C ) we see, that due to £4 = two of the terms in Bi(t) combine 
together; this corresponds to the fact that T-b 3 (lvi) has only 7 weights while Td 4 (wi) 
has 8. Analo gous, but more complicated combinations and cancellations take place in 
the proof of Wlh. 
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Example 3 The case with G 2 can be obtained from D4 after imposing a symmetry 
condition with respect to the outer automorphism Vi (3.41) of D 4 of order 3. The 
restrictions that V\ imposes on the scattering data are: 

Ci = 0, Cl-C2 = C3, <A>4 = 0, (p 01 - lfi 02 = (fi03- (3-55) 

or in other words 

Vl(0 = C, «l(#)) = $D- 

Then one can check that due to ( p-4^j ) i\t) *s invariant with respect to v\: Vi(q(t)) = 
q(t) and, consequently, is an element of the subalgebra G 2 . 

Indeed, if we average the set of simple roots 0/D4 with respect to the action of ' v± 
and also the system of equations (3.4< ) then we get the system of simple roots of G 2 : 

Pi = ^(ai + a 3 + a 4 ) = -(ei - e 2 + 2e 3 ), [3 2 = a 2 = e 2 - e 3 , 
and the system of equations 



(3.56) 



Qi.tt 



-(Q,/3i) 



(3.57) 



which is slightly different from (1.1) for G 2 ; namely the variables, let say q[, q' 2 and 
q' 3 = —q[ — q' 2 , entering in the original system (1.1) for G 2 are 

q[ (t) = Q t (t) - 2 In 3 q' 2 = Q 2 - In 3 (3.58) 



The solution to (3.57) is provided by: 
Qi = lnBi(t;C, <pb), 



Q 2 = In 



B 2 (t;(,<p ) 



(3.59) 



Bi(*;C,^) 

where Bi^ 2 (t) are obtained from ( 3.44 c ), ( 3.44b ) with the additional restrictions ^3.5$ ); 
namely we have: 



B l (t;t,vh) = 2b b 1 

B 2 (t;(,Vo) = b b 2 1 
, 3(Ci +C2) 



■ cosh<^i 



C1+C3 , C1+C2 , 

■ cosh tp 2 H cosh 953 



Ci C2 ~ r ~ ' Cs 

d 2 cosh(0i,^) C| cosh(31,<#) C| cosh(32,<£) 



26?, 



C2C3 C2 — C3 



C1C3 Ci + C3 C1C2 Ci + C2 



C1C2 



Ufin -*\ 1 3(Cl+C3) U f--i . 3(C 2 - C3) ,,„-i 

cosh(10, y>) H — — cosh(ll, tp) H — — cosh(21, <p) 



C1C3 

5o 



8(Ci + C2)(Ci+C3)(C 2 -C3)' 



C2C3 



61 



(3.60) 



8CiC2C 3 ' 



whe re by ij we denote the root i(3\ +jf3 2 ofG 2 . Obviously each of the terms in Bk(t) 
in ( 3.60J can be related to a weight of the corresponding fundamental representation 
T(u>k) of G 2 ; so again we cast the solution in the form ( 1.3<\ ), (1.31). 
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We summarize the results of this section by the following remark. One can view 
the solutions of the CTC related to the classical series B r , C r as obtained from the 
Moser formulae combined with the corresponding constraint on the scattering data. 
The case of D r can be obtained likewise if we take the Lax matrix to be pentadiagonal 
as in (|E§). 

On the other hand starting from a Lax matrix related to each of this series we can 
always apply Moser's approach and derive as a result the functions qk(t), k = 1, . . . , N. 
In order that both answers for the CTC-solutions be compatible one needs to show 
that 

Qk(t) = -q~ k (t), k = N + l-k (3.61) 
where N is the dimension of the typical representation of g. Equation ( |3.61 ) can be 



derived from the results in Sections 2 and 3 and shows the compatibility of the two 
approaches to the CTC solutions. 



4. Dynamical regimes and large time asymptotics 

There are important differences between the RTC and CTC, especially the asymptotic 
behavior of their solutions. Indeed, for the RTC, one has 0, [l2) that, both the 
eigenvalues, and the constants </?/c(0), are always real- valued. Moreover, one can 
prove that (j, ^= Q for k j, i.e. no two eigenvalues can be exactly the same. As a 
direct consequence of this, it follows that the only possible asymptotic behavior in the 
RTC is an asymptotically separating, free motion of the particles. 

The situation is different for the CTC. Now the eigenvalues Cfc = K fc + i-Vk, as 
well as the constants fk{0) become complex. Furthermore, the argument of Moser 
[Q does not apply to the complex case, so one can have multiple eigenvalues. The 
collection of eigenvalues, (k, still determines the asymptotic behavior of the solutions. 
In particular, it is Kk that determines the asymptotic velocity of the k-th particle. For 
simplicity, we assume ^ Q for k ^= j. However, this condition does not necessarily 
mean that Kk =/= Kj. We also assume that the k^'s are ordered as: 

Kl < K2 < • • • < UN- (4-1) 

This ordering is known as the sorting condition. More generally it can be understood 
as — k e Wd - the closure of the dominant Weyl chamber. Once this is done, for the 
corresponding set of ./V particles there are three possible general configurations: 

A) free particle propagation (Moser case) - then qk (t) have linear in t asymptotic 
behavior and — k is in the interior of Wd- 

B) bound state(s) and mixed regimes when one (or several) group(s) of particles 
form a bound state - then each group of particles oscillate around common trajectory 
with linear in t asymptotic behavior; then — (k, a^) = for some set of indices k S lbs- 

C) degenerate solutions when two (or more) of the eigenvalues = (k+i = ■ ■ ■ 
are equal - then quit) — qk+i(t) have logarithmic in t asymptotic behavior, i.e. the 
distance between the particles grows as Int. 
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Obviously the cases B) and C) have no analogues in the RTC and physically are 
qualitatively different from A). 

4-.1. Asymptotically free regimes 

We begin with the first possibility from the above mentioned - free particle asymptotics 
(Moser case) . It is realized if we require that all real parts of the eigenvalues Kk are 
pairwise different; i. e. — k belongs to the interior of the dominant Weyl chamber Wd- 

(-£,a s )> 0, s=l,...,r, (4.2) 

while the imaginary parts may be arbitrary. 

Let us now consider the asymptotics for the A r , B r , C r and D r series. Using the 
explicit expressions for Bfc(i) it is not difficult to evaluate their asymptotic behavior 
for t — > ±oo: 

Etas® = ^'(^le'-**^ (1 + 0(eMTK±t))) . (4.3) 

Here are the highest (lowest) weights, related through the Weyl group element wq: 
uifr = wq(uj^) and 

K+ = min Re (<f, 7 - ) = _ («> a k), 
7 erv+)V+ 

= min Ra(C > 7-Wfc) = -(/t,a i ) 1 (4.4) 
7 er( w +)\a ) - 

see Appendix B. 

Note the natural way in which the two asymptotics are related by the wq 
transformation of the Weyl group, namely: 

= «*> (B+Jt)) = W( fc )(C>oK + ))e ( ^ t) ^K + )) 

= W^(wo(C),Uk) e{wam)) ' ui) - ( 4 - 5 ) 
More generally since W^ k '{Q, 7) depends only on the scalar products of the type 
(O7) the action of the automorphism wo on W^(Cil) is given by (£, w>o(7)) = 
W^(w ((),j)- 



The relations (4.3) and ( |4.4| ) are due to the simple fact that the leading exponent 
for t — > 00 (t — > —00) in Bk(t) corresponds to the weight 7 £ T(ujk) for which the 
value of Re (— £, 7) is maximal (minimal). Since — k € Wd this maximum (minimum) 
is realized when 7 = ui^ (respectively, 7 = ui^ ) . 

From the previous considerations we have: 

Qk{t)=ln- ? - s = \ - v In B s {t). 4.6 

Bk-i(t) ^-J {a s ,a s ) 
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and consequently the asymptotics g~ of q(t) for t — > ±00 are given by 

£(t) = -20E + + £ T^H InWW(C, (4.7) 

= «*> (-205 + 9*>) + E t-^tt In ( Wo(C ), w +), (4.8) 

up to terms falling off exponentially for t — > ±00. The explicit expressions for the 
components qk(t) for the A r series are well known, see e.g., M O, ^]. For the other 
classical series of Lie algebras we obtain: 

q^ m (t)=TKkt±^0k+Pk + o(e^ t y k = l,...,r, (4.9) 

with 

h = In (w k TT (2C S - Cfc) 2 ) , ^± = min Kf. (4.10) 
The only exception to ( [4.9|) , ( |4.1C ) is for g ~ D r with odd r, k = r and i — > —00; 

then 

?na.(*) = ~ 2 Cri + <y50r + In ^tvr n( 2 Cr + 2Cs) 2 ^ . (4.11) 

It has been known for a long time 15, [l6| that for the RTC the asymptotic 
velocities iP 1 are related by v~ — wo(v + ). In the complex case the analogs of v are 
the complex vectors — 2£ and — 2u>o(C) respectively. 

Up to now we know of only one physical application of CTC as a model describing 
the iV-soliton train interactions |l7], [l8|, [l9|, ^0). Getting insight from it we will 
interprete Keqk(t) as the trajectory of the center of mass of the fc-th 'particle' (soliton). 
Besides each particle is complex and possesses an internal degree of freedom. Then 
— 2Rc(fc = — 2k/j will be the asymptotic velocity of the &;-th particle at t — > 00, while 
— 2ImO = — 2r]k determines its asymptotic phase velocity. 



4-2. Mixed regimes for B r , C r and D2 



Our aim in this and next subsection will be to consider the cases when two or 
more particles form bound state(s); we will say that several particles form a bound 
state if they have equal asymptotic velocities. In this subsection we consider only 
those members of the classical series for which wq = —id. 

Bound state(s) are possible when — k is on the boundaries of Wd, i-e. if we have 

-(a k ,K) = 0, fee lbs, (4.12) 

where lbs C {l,...,r} is a subset of indices. If I^s = { m } contains just one 
index m < r then n m — k,„+i and we will have a two particle bound state. If 
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-^bs = { m , m + 1, . . . .m + p} and m + p < r then n m — K m +i = . . . = n m+p and we 
have ap + 1-particle bound state. The cases when the largest index in ibs is equal to 
r should be considered separately; indeed due to the fact that the sets of simple roots 
for the different classical series differ only in the choices for a r this cases may lead to 
substantially different results. 

In our previous paper pp| we obtained the large time asymptotics for the two- 
particle bound states in the A r CTC. Here we will briefly analyze more general cases 
when: a) g belongs to the other classical series and one bound state may be present, 
i.e. when lbs = { m }i m < r and lbs = {m, m + 1}, m + 1 < r; b) two bound states 
may be present, i.e. ibs — { m iP}i m + I < p < r. For brevity we will write down 
the asymptotics only for those components qk(t) which differ from the typical ones 
(4.9). We will limit ourselves with the cases when the mixed regime contains two- and 
three-particle bound states only. The other more complicated regimes can be analyzed 
analogously. 

Indeed, if for k G lbs we have (k, oik) = then at least two terms in Bk(t) may 
have the same asymptotic behavior. To our purpose it is sufficient to evaluate only 
the leading exponents. Thus we find 



B± (*) = e 



M«),<) 



where 

K'* = 

p 



mm 

7 er p , ± o+) 



E 



W {p \c 



O 



(ft**) 



(4.13) 



-<*><*).«) W(p)(£ w ± t a), 



Gp(fi) = < a > 0, (a, k) 



0, ± )' pl > 1 
(a, a) 



(4.14) 



r p ,±(^+) = T(iv+)\ {wp ,Wp Ta,a e G± (re)} . 

The condition (a, k) = ensures that only oscillates when t — > ±oo while 

the third condition in Gp (k) means that lo^ =p a € r(w+). 

We start with the simplest case lbs = which contains several qualitatively 

different subcases which will be listed below. In each of them it is possible to describe 
the sets of roots G^(k) and to evaluate the estimating exponents K'^ 1 , for details 
see Appendix B. It turns out that G^{k) = for p ^ m and G+(k) = {a m }, 
G~(k) = G+(u>o(«)) = {wo(a m )}. In what follows we will concentrate mainly on 
the asymptotics for t — > oo; the asymptotics for t — > — oo we will obtain by formula 
(4.8). Therefore the asymptotics of B^, AS (t) for p ^ m will be given by (4.3) while for 

#m,as(*) we § et: 



s+ (*) 



-.-K'+t 



(4.15) 



W^ m \C,u+) + W^\a m )+0(e 

with the following result for K'^ valid for any of the algebras in the classical series 
(see Appendix B): 

K'+= min [-2(K,a s )\. (4.16) 

s:(a s ,Q m )/0 
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In other word the minimum should be taken with respect to the simple roots ot s that 
are connected to a m in the Dynkin diagram. 

Now we are in position to compare the asymptotic velocities of the particles and 
to single out the structure of the bound states (if any). Since the asymptotic velocity of 
quit) for t — > oo is equal to — 2K& we just have to see what constraints on {ki, . . . , K r } 
will be imposed by — (k, at) > 0, k ^ m and — (k, a m ) = 0. For g ~ B r , C r and 
m < r we have 

Ki < ... < K m = K m+ i < . . . < K r < 0; (4-17) 
for g ~ T>2n and m < In — 1, to = In — 1 we have 

«1 < . . . < K TO = K m+ i < . . . < K 2n -1 < — |«2nl < 0, (4.18) 

and 

Kl < . . . < K m < K m+ l < . . . < K 2n -1 = K2n < 0, (4-19) 

respectively. 

Finally for g ~ B r , C r and m = r we get: 

Ki<...< K r -i < k,. = 0, (4.20) 

and for q ~ D2« and m = 2n 

«!<...< «2n-l = ~«2n < 0. (4-21) 

From ( 4.17j )-(4.1S) it is easy to see that for m < r we always have one bound state 
of two particles (m-th and m+ 1-st); the rest of the particles go into a free asymptotic 
regime. If in addition wo — —id , as we assumed in the beginning of this subsection, 
this bound state will be present also for t — * ~ oo. Therefore for this classes of algebras 
we have stable 2-particle bound states for all m < r. 



For m = r the situation is different. From ( 4.20| ) we see that the condition 



— (k, a r ) — for g ~ B r , C r just means that the r-th particle has vanishing velocity. 
As for g ~ T>2n the condition — (k, a.2n) — means that the 2n — 1-st and the 2n-th 
particle have opposite velocities. Therefore for m = 2n no bound states are possible. 

The next possibility is that the set lbs = There are qualitatively different 

cases here: a) (a m , a p ) = and b) (a m , a p ) ^ 0. Each of the values m and p in case 
a) can be considered independently and to each of them applies the analysis already 
posed above. In the generic case m < p < r we will have two pairs of bound states 
each containing two particles; if m < p = r then we have only one bound state of two 
particles. An exception here is the case g ~ T>2n and m = 2n — 1, p = 2n. The two 



roots OL<i n — i and ct2n are obviously orthogonal, but now the condition (4.12) leads to 
K2n = K-2n-i = and as s result in this case we have only one bound state consisting 
of two particles with vanishing velocities. 

Let us now analyze the case b). For generic values of m < p < r what we find is a 
bound state of three particles. One possible realization of b) is to take p = m + 1 < r; 
then the condition ( 1.12| ) leads to 



«!<...< K m — Km+l — K m+ 2 < . . . K m +3 < ■ ■ ■ , (4.22) 
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i.e., the particles numbered by m, m + 1 and m + 2 move with the same asymptotic 
velocities and form a bound state. Again we must look through all possibilities when 
case b) takes place and point out possible exceptions. Such for example is the case 



when q ~ B r , C r and m = r — 1, p = r. Equation (4.12) then gives K r -i = K r = 0, 
which means that this is bound state of two particles: r — 1-st and the r-th with 
vanishing velocity. 

If g ~ D2n and m = 2n — 2, p = 2n — 1 we get a three-particle bound state 
with velocity K2 n —2 — K>2n-i — K 2n < 0. The last example related to this algebra is 
m = 2n — 2, p = 2n which corresponds to /t2n-2 = ^2«-i = — K 2n- This means that 
the particles 2n — 2 and 2n — 1 form a bound state, but the last 2n-th particle moves 
with the opposite velocity and is not a part of the bound state. 

Obviously the number of examples can be extended to include sets lbs with more 
indices; one can expect to have bound states with increasing number of bounded 
particles. It is not difficult to present also the explicit form of the asymptotics 
of qtasify- T ne mos t difficult part in this calculation is to determine the sets of 
roots Gp (k). We list these sets of roots in Appendix B for the classical series of Lie 
algebras related to the sets lbs with one and two indices. Indeed if choose g ~ B r , 
-fbs = {r — 1>7"}- Then the sets of roots 

GJ-\{k) = { a r- 1) £* r _i + a r , a r -i + 2a,.}, Gf (k) — {ot r , ov-i + a r ,Q r -i + 2a r }, 
Then 

+ W^-V (cv-i + Or) + W^-V (a r _i + 2ar)l , (4.23) 
B+ &s = e^W'^) \w ir \(,Lu+) + W {r \a r ) + W {r \a r ^ + a r ) + W (r \a r ^ + 2a r ) . 



Now we have to insert ( 4.23 ) into (f4.6| ). After some calculations we obtain 
K r = and 



= "2i*-l* + + In W( J;^ _ 2) 

g+ as (<) = -2 l?7r < + <p 0r + In — (4.24) 

S P = WW(C,0 + S P , p = r-l,r, 

S p = W {p) {a p ) + Vy (p) («r-i + otr) + W {p) {a r ^ + 2a r ). 

In this subsection we presented various types of mixed regimes which could be 
called regular. By regular here we mean that the number and the structure of the 
bound states at t — > — oo coincides with the ones for t — > oo. In the next subsection 
we consider the 'irregular' mixed regimes, which change qualitatively their structure 
during the evolution. 
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The 'irregular' mixed regimes take place only for algebras for which wq 7^ —id. This 
takes place for g ~ A r and g ~ D271+1. At the end of this subsection we will explain 
why the 'irregular' regimes, i.e., the effects of 'creation' and 'decay' of bound states 
can be related only to q ~ T>2n+i- 

First of all we note that most of the bound states related to T>2n+i are regular. 
Such are the states corresponding to lbs = {m} with m < 2n and ibs = {m,p} with 
m < p < In. The formulas for the asymptotics in all these cases are quite analogous 
to the ones already presented. 

Let us start with the first 'irregular' case with = {2ri}. This means that 

K 2n = K 2n+1 8£Qd f°r t ~ * ioo 



92n.as 



T2K 2n t =F 2ir/2„t ± <^ ,2n + ^'(C)> 

itn+l^s = -^2nt - 2ir/2n+lt + ¥>0,2n+l + #to+l(0> ( 4 - 25 ) 

M 2 ")(C,<) + W^(a 2n )] W^ 2 " +1 >(C,< +1 ) 
/3 2 + /(C)=ln i 

/?2 + ^i(C)=ln 



^ (2 " +1) (C,^2 + n + l) 

W^){lut n ) + W^){a2n) 



and /3Z''(C)i k = 2n, 2n + 1 are obtained from P^''(C) by using (1.5). Obviously at 
t — > —00 the 2n-th and the 2n + l-st particles have opposite velocities, while for i — ► 00 
their velocities become equal. This situation can be viewed as 'creation' of a bound 
state. 

The second 'irregular' case is with lbs = {2n+l}. This means that K 2n — —K 2n +i 
and for t — » ±00 



^2n.as 



Pt"(C), 



= T2kW T 2l?72n* ± ¥>0,2n 
qfn+l.as = 2k 2«< - 2il]2n+lt + </?0 : 2n+l + Ain+1 (0) 

^ 2n+1) (c> 2 Vi)+^ (2n+1) K + i; 



/?2 + /(0=ln 



/^+l(C)=ln 



(4.26) 



^(2»-D (C,<_r) 

^ 2 " +1 KC,< + i) + ^ (2n+1) (« 2n+ i) 



and again /3I'\Q, = 2n, 2rt + 1 are obtained from /3t''(Q by using (1.5). Now at 
i — > — oo the 2n-th and the 2n + 1-st particles have equal velocities, while for t — > oo 
their velocities become opposite. This situation can be viewed as 'decay' of a bound 
state. 

The next more complex situation is when Zbs = { m ,p}- Again we should consider 
two distinct subcases, namely (a p , a m ) = and (a m , a p ) ^ 0. 
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In both cases we recover 'regular' asymptotics provided m < p < 2n; namely 
for such choices of I\> s we have either two pairs of two-particle bound states (if 
(a m , a.p) — 0) or a three-particle bound state (if (a m , a p ) ^ 0). 

The 'irregular' cases with (a m , a p ) — are two types. The first one takes place if 
m < 2n and p = 2n (p = 2n+l). Then for t — > — oo (t — ► oo) we have two bound states 
formed by the particles {to, m + 1} and {2n, 2n + 1} while at t — > oo (i — > — oo) the 
second bound state decays and we are left with only one bound state. Quite different 
is the situation when — {2n, 2n+ 1}. This corresponds to n 2n = n 2n +\ = 0, so this 
is a regular case but with only one bound state formed by the particles {2n, 2n + 1} 
with vanishing velocity. 

There are only two 'irregular' cases with (a m , a p ) ^ 0, namely I^s — {2« — 1, 2n} 
and lbs = {2n — 1, 2n + 1}. The first one leads to K2 n —i = K 2n = ^2n+i < and to 
the following asymptotic behavior of quit), k = 2n — 1, 2n and 2n + 1: 

^-i,aB = T2K 2n -it T 2in 2n ~it ± <A),2n-i + /^'n-D 

92n,as = T2K2n-l£ T 2i?72„i ± </? ,2n + /^'n > ( 4 ' 27 ) 
?4+l,as = -2K2n-l* - 2in 2n +lt + V0,2n+1 + /3 2 'n+l> 

) ( c,<_ 2 )' /32 "~ ln e 2 .„+_ i( c) 

_ ln ^ (2 " +1) (C,< +1 ) 

^2n IU 

S;-+(C) = W W (C,w+) + W^(a p ) + W^(a 2n ^ + a 2n ). 

From these formulae we find out that for t — > — oo we have a two particle bound 
state formed by 2n — 1-st and 2n-th particles, while for t — > oo the 2n + 1-st particle 
'joins' them and we have a three-particle bound state. 

The case with I^s — {2n — 1, 2n + 1} is analogous: the only difference is that 
at t — > — oo we have a three-particle bound state formed by the 2n — 1-st, 2n-th and 
2n + 1-st particles, while for t — > oo the 2n + 1-st particle 'separates' from them and 
we are left with a two-particle bound state. 

Let us analyze this situation on the bases of our remark at the end of Section 3. 
Let us first explain why such irregular solution are not possible for the A r -series. In 
this case we have r + 1 particles and the sets of asymptotic velocities for t — > oo 
and t — > — oo differ only in the ordering: {— 2K r+ i < —2k,. < . . . < — 2ki} and 
{— 2ki > —2k 2 > . . . > — 2k,. + i} respectively. That is why it was quite natural to 
identify the fc-th particle at t — ► — oo with the fc-th particle at t — ► oo: they move with 
equal velocities. This is compatible with the action of wo in the A r case, see ( |B.2| ). As a 
result if we have, say two bound states at t — ► — oo, i.e. — 2k± = —2n 2 > — 2«3 = — 2«4 
at t — > oo we will have again two bound states — 2/t r +i = ~2n r < — 2/t r — l = — 2ft r _2- 

Next we can view the solutions of CTC related to the clas sical series B r , C r and 
D r as special symmetric solutions of the sl(N)-CTC, see (3.61). Then it is enough to 
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consider only 'half of the trajectories; the other half is obtained as a 'mirror' image. In 
this situation the sets of initial and final velocities are different and the identification, 
good for A r are not possible for the other classical series; quite different is also the 



action of wq on the orthonormal basis of the root space, see (B.2). 

As we mentioned above, if we consider the whole picture with all the N- 
trajectories we will see that no 'creation' or 'decay' of bound states take place. In 
the cases of D2n+i and (k, a^n) = (or (k, ctan+i) = 0) this can be explained as 
follows. For t — > oo (t — * — oo) we have a bound state between the 2n-th and 2n + 1-st 
particle which for t — > — oo (t — > oo ) transfers into a bound state between the 2n-th 
and its 'mirror' symmetric In + 2-th particle. 'Cutting' off the symmetric trajectories 
with numbers 2n + 2, . . . , N = An + 2 we find the effects described above. 

Analogously we can explain the situation with a three-particle bound state at 
t —* ±oo and a two-particle bound state at t — > =Foo. The whole picture of An + 2 
particles will always contain a three-particle bound state. 

4-4- Bound state regimes. Periodic and singular regimes. 

These regimes take place if k = 0, i.e. the set of eigenvalues Cfc = ink are purely 
imaginary. Then each of the functions Bk(t) will be generically bounded. In particular 
this means that all the complex 'particles' (or solitons) will move together forming a 
bound state with a large number of degrees of freedom. In order to avoid degeneracies 
we have to request that fj £ Wa- 
in order to have periodic solutions we need one more restrictions upon r/k , namely 

Tjk ~ Vm = Sfe m r? , (4.28) 

where Skm are integers; if Skm are rational we can always make them integers by 
rescaling 770- 

Example 4 Let g ~ B2 = so(5). The corresponding equations have the form: 

-*»-», ^| = - e *-9i +e -*, (4.29) 



dt 2 ' dt 2 

and their periodic solutions are given by: 



B 2 (t\ 

gi = lnflift), 92 = ln-^, ( 4 -30) 



where 



n M 1 / cos2 Pl (ci>(t)+r) cos2 P2 (<i>(t)-r) p\-p\ , 

6iW = i6^?-p!)1 v\ + Pi + (431a) 



m = fcos^w+p-r) + cos( P ^(t)+P + ry } 

&V0P1P2 { p+ p- 

where rjk = PkVo, Pk o,re integers and 

$i =7/ i+- + , r = - , p±= Pl ±p 2 . 4.32) 

4 V Pi Pi J 4 V pi p 2 
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The period is provided by 

t = — — — (4.33) 

where so is the greatest common divisor of p\, p2, p+ and p- . 



Our next remark is that in the generic case when Re0ofc 7^ the solution (4.31a 



(4.31?) is a regular one; then |£?i(i)| and |£?2(i)| are strictly positive for all t. If we 
choose however, Re0oi = Re0o2 = then |$i(£)| and |i?2(i)| may vanish and the 
corresponding solutions qi(t) and 92 (t) become singular. Due to the periodicity, if 
|£>i(t)| and |£>2(t)| vanish at certain points tgi and io2 respectively, then they will 
vanish also at toi + kr and £02 + kr for any integer k — 0, ±1, ±2,. . . . 

Example 5 Let g ~ C2 — sp(A). Since the algebras B2 — C2 then the corresponding 
solutions differ by change of variables. Let us denote all variables of the C2-CTC 
model by the same letters as for B2, adding additional 'bar' to distinguish between 
them. Then the C2 -CTC system has the form: 

4^ = e" , 4f = -e M + 2e- 2 ' 2 , (4.34) 
at at 2 - 

and the solution is presented by: 

gi(t)=]nBi(t)+ln2 J q 2 (t) = In + i In 2, (4.35) 

where 

B 1 (t) = 2B 3 (t,Ci.C2,0o 1 ,^oa), = 2B x (t, fi, C2, 0oi, 0o 2 ), (4-36) 
and Bk(t) are given by( J h .31c ), ( 4-31b\ ) and 



Ci — Ci + C2, C2 = Ci - C2, 

001=001 + 002 , 002 = 001 - 002- (4.37) 

0/ course these last two examples are analytic continuations of the solutions 
presented in J71|/. 

In analogy to the previous example we may assume ^ to be purely imaginary 
with r]k = pkijo with integer p^- Then we get the corresponding periodic solutions 



to Eq. (4.34). More generally, inserting purely imaginary values for will result in 
periodic solution in all the above examples. These solutions may become singular if 
the corresponding parameters 0ofe are purely imaginary. 
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4-5. Degenerate solutions 

Let us briefly discuss the degenerate solutions to the CTC. The degeneracy is possible 
only if the matrix 1/(0) has nontrivial Jordan cells p0| . 

One possibility to derive the de gene rate solutions is to evaluate the limit £i — ► 
(2 — * • • • — ► (k of the solution (3.9), ( 3.10 ) using the l'Hospital rule. 

If in particular, we have complete degeneracy (i.e., all £fc are equal to zero) the 
solution of the sl(N)-CTC can be obtained in a more simple way. Since all -Bfe(i) are 
polynomials in t which must satisfy 

B (t) = B N (t) = 1, B k B k -B\ = B k ^B k+1 . (4.38) 

we find that they depend on N — 1 constants fk, k = 1, . . . , N — 1. More specifically 
B k (t) must be polynomial of degree k(N — k) whose coefficients can be determined 
explicitly for example, by the method of undefined constants. 
For example, for TV = 3 and Ci = C2 = C3 = we get: 

Bi{t) = -\f + fit + /a, B 2 {t) = -\t 2 + fit - ft - /a, (4.39) 

and -B3 = 1, where f k , k = 1,2 are complex constants. If fi is real we may use the 
translational invariance of the CTC-equation and change t — ► t + /1 to eliminate it; 
then the solution ( 1.39| ) can be written in the form: 

S 1 (t) = -i* 2 + Fi, B 2 (t) = ~t 2 -F u Fl = l -f 2 1+ h. (4.40) 



If F\ is real, then the solutions q k (t) have singularities for t = ± \J 2 1 F\ \ . The large 
time asymptotics are given by: 

lLs(t) = ~qf, as (t) = 21nt - In 2 + ivr, ^(i) = 0, (4.41) 
i.e., they do not depend on the constants f k and are complex. Analogously for N = 4 
and £1 = C2 = C3 = C4 = we find: 

B (t) = B 4 (t) = 1, B,{t) = \{t 3 + ht 2 + f 2 t + / 3 ), 

6 

B 2 (t) = - l -ht* - i/ft 2 - i(/ x / a - 3/ 3 )t - - 2/i/a), (4.42) 

B 3 (t) = - i/it 2 - {\ft - \h)t - ^(2/f - 6A/ a + 9/3), 

where A; = 1, 2, 3 are complex constants. If f\ is r eal w e may change t — > t+fi/3 to 
eliminate one of these constants; then the solution ( 4.42| ) can be written in the form: 

Bi(t) = \(t 3 + Fit + F 3 ) , B 2 (t) = ~ (M 4 - 6F 3 t + F 2 ) , 

B 3 (t) = ~ (t 3 - F 2 t + F 3 ) , (4.43) 
where Fi and F 2 are expressed through f k by 

F 2 = /a - ^ , F 3 = / 3 + ^ - M. (4.44) 
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Obviously, these solutions will be regular if B k (t) have complex roots and will 
develop singularities if one (or more) of their roots are real. More specifically if F 2 = 
the solution becomes symmetric, i.e. B\(t) — —B^{t) and has a singularity at t = 0. 
If in addition F3 is real, then there are singularities also for t = — and t — 

The asymptotics of these solutions are easy to calculate: 

qi, as ±(t) = 31nt-ln6, «a,««±(t) = Int - In 2 - m, (4.45) 

33 >M ±(*)=-In* + ln2 1 q4,as±(t) = -31nf + hi6 + i7r. (4.46) 

Note that again these asymptotics: a) do not depend on the constants F k , and b) are 
always complex. The last property is a consequence of the fact that degeneracy is 
possible only for the CTC. 

5. Conclusions 

Detailed analysis of the properties of the fundamental representations of the simple 
Lie algebras allowed us to propose an effective and invariant parametrization for the 
solutions of the CTC. These solutions describe much richer asymptotical regimes as 
compared to the RTC. The explicit solutions proposed above allow one to evaluate 
explicitly the large time asymptotics for the whole variety of dynamical regimes. 
Further investigation deserve also the degenerate solutions. 

One final remark is that one more step is necessary for the perfection of the 
) , namely one should look for an invariant expression for 
the functions ((f, 7/*^) hi terms of 7^ and the roots system A only. Work in this 
direction is in progress. 
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explicit formulae (Oq), (1.36 



Appendix A. The properties of V. 



Here we outline some of the details in deriving the expressions for A k (t) and Bk(t). 
As we mentioned in Section 2 we need the explicit expressions for the minors of V. 

Let us consider the eigenvalue problem (^) and let us make use of the explicit 
tridiagonal form of L(0). Then it is not difficult to find that the eigenvector related 
to Cfe is of the form: 

/ r k p x \ 

r k P2{(k + Po(Cfc)) 

M 



.,(fc) 



rkPsia+PiiCk)) 



(A.l) 
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where P s (Ck) stands for a polynomial of degree s in ( k and 
1 1 



Pi = 1, 



P2 



P3 



Pn 



(A.2) 



ai(0)' " ai(0)a 2 (0)' ~" oi(0) . . . ajv-i(O) 
Next we note that the terms with P s (Ck) do not contribute to the minors 

V 



where 



k 



and so 
V 



k 
ik 



■ n k pi . ..p k W(ii, . . .,i k ), 



,4fcJ 



det 



1 

2Cn 



1 



\k-l 



n -c P ), 



(A.3) 



(A.4) 



(2C il ) fc - 1 ... (2Q hJ 

is the Vandermonde determinant and / = {ii, . . . , ik}- Next we have to take care of the 
factors p s , which can be expressed through <f(0) since afe(O) = | exp(— (q(Q), afe)/2). 
We remind also that are determined up to a sign by the normalization condition 
(2.4). These remarks and the properties of the fundamental representations of the 
series A r (2.19), ( |2.2C ) are sufficient to treat the A r series. 

Let us now derive the symmetry relations (2.9), ( [2.10 ) for the B r and C r algebras. 
To this end we introduce the matrices S as follows: 



S = E(- X ) fe+1 ^k + E kk) + (~l) r E r+1 . r+1 
k=l 

r 

= J2(-l) k+1 (E k - k -E- kk ) for C r , 



for B r , 



k=l 



E(-!) fe+1 (^ 



E kk) 



(A.5) 



for D r , 



k=l 



which enter into the definition of the corresponding orthogonal and symplectic 
algebras. By Ej k we denote a n N x N matrix whose matrix elements are equal 



to (E jk ) 



jk Jmn 



5j m 5kn and as in ( |2.9| ) k 



N + 1 — k. Then we make use of the fact 



that if V is a group element of the corresponding group then V T 

Tk =V ■ 



SV^S- 1 , i.e. 



1 . . 
1 . . 



k 
k 



N 
N - 1 



(A.6) 



where the 'hat' means that the index k is missing. Equation (A.3) and (A.6) readily 
give: 

W(l,...,k,...,N) A 



r kT k 



p- k W(l,...,N) 



n 



P s r s 



(A.7) 
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Taking the product of (A/?) for k = 1 to r one is able to evaluate n«=i r k m terms of 
p s and Q s alone. Next putting p s , N and £ s as appropriate for the series B r and C r 
we derive the expressions for w k ( 2.11 ) , ( 2.12 ) and ( 2.13 ) . 
For the D r series the matrix L(0) is of the form: 



ai b 2 



L(0) = 



V 

For the fc-th eigenvector we get: 



b r —i a r ~i 

CL r —l b r 



a r 






ft,. 



a r 






ft,. 



— b r ftr — 1 

ft r _i — br—1 



-b 2 a x 
oi -61/ 



(A.8) 



r kPr(C k 



/-r~2 



- C/c fc ) 

rkQPr+2 



(A.9) 



where C is a coefficient to be calculated below. The symbol ~ in ( A.8 ) means that in 
the right hand side we have omitted terms polynomial in which do not contribute to 
the minors of V. Note also that C enters into play only when we need a minor of order 
r or higher. Such necessity appears in two cases: when we evaluate the expression for 
A r (t) and when we derive the symmetry relation. 



From (2.3) and the explicit formula for L(0) (A.£) we find that p^ in ( A.9 ) are 
given as follows: 

fc-i 



Pi 



1. 



Pk 



n 



1 



Pr = 



1 



Pr+l = 



n 



for k 



1, 



2a r (0) a s (0) 



Pr+2 = 



1 - 



\ o.(0) 



Pr+fe = Pr+2 



fc-1 

n 

s=2 



,(0)' 



for fc = 3, 



(A.10) 
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The determinant of V gives: 

2r 

1 = detV = 2(-l) r - 1 CW{l, . . . ,2r) TT (A.ll) 

„_1 Cs 



In analogy with flA.q ), (A. 7) we get: 



2c 



(-l) r 2C<; p W(l,...,p,...,2r) fjp s r 



p 2r W(l,...,2r) 11 C 



n^p. (A.i2) 



Taking again the product Jl P =i r v r v m and substituting the expressions for p s 

from (A. 10) we find 

r 

c=( _ 1)r+ l"Q Cs) (A13) 
s=l 

and in addition - the relation ( [2.14| ). 

Now it is easy to find the expression for the minor of order r: 

V { I ■ ■ I } = \ (l + 1^1) . • • , ir) jft n s P, (A.14) 

needed for the derivation of A r (t) ( 3.32| ). 
Appendix B. Algebraic details 

The action of Wo on the simple roots is well known |m| : 

w (a k ) = -a~ k , (B.l) 

where k — r — k + 1 for A r ; k = k, k = 1, . . . ,r for B r , C r and D2„. For g ~ D2 ra +i 
we have fc = k for A; < 2n — 1, and u^c^n) — — o^n+i, wo(a2n+i) = — o^n- More 
specifically wq acts on the orthonormal basis {e k } in the root space as follows: 

w (e k ) = e k , forA r , (B.2) 

wo(efe) = -e k , for B r ,C r ,D 2 „, 

and for T> 2n +i 

w (e k ) = -e kl for k = l,...,2n, wo(e2n+i) = e 2 „+i- 

Next well known fact is that the weight system r(w) is determined uniquely by 
the highest weight u>. The reconstruction of the weights 7 S r(w) is performed by 
using two facts: i) if 7 G r(w) then 1^(7) G where w is any element of the Weyl 
group; besides 7 and 10(7) have equal multiplicities; ii) if a > is a positive root and 

^4^>0, (B.3) 
(a, a) 

then to — sa G r(a>) for all s = 1, . . . ,p. 
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In particular, if u = u>t and a = Y2 s =i m s a s we see that (B.3) is fulfilled only if 
m k — p > 0. Thus we find that generically (i.e. for k < r) along with ut, weights in 
r(a-i^) are also 

72 = u+ - (ctk-i + a k), etc. 



7i =^ k 



Oik-, 



Using the sorting condition (4.2) we easily find that 



mm 



7 er«)W 



[-(«,^ + -7)] = 



-(re,a fe ),(B.4) 



which proves the estimations in Equations (4.3). 

The same method can be applied also when one of the roots satisfies (re, a m ) = 0. 
As a consequence of this condition at least two terms in Bk(i) may have the same 
asymptotic behavior. 



Here we will first describe the sets of roots G+(re) see (4.14) and then will also 
outline the proof of ( 4.16| ). Obviously if (re, a m ) = only G+ (re) will be nonempty and 
will coincide with {a m } and therefore r Pj+ (w+)\{a;+} while r mi+ (w+)\{a;+ , w+ — a} 
where a — a m + ^2 s m s a s . The minimum of — 2(re, w+ — a) will be achieved if we 
limit ourselves with roots a of height 2. Now it remains to take into account that 
a m + a s is a root if and only if (a m ,a s ) < 0. The corresponding result for t — > — oo 



is obtained by acting with wq. This proves (4.16) 



We finish this appendix by describing the sets of roots Gp(re) for (re, a m ) = 
(re, a m ) = 0. First, if (a m ,a p ) = then G+(re) = {a m }, G+ (re) = {a p } and all the 
others G^(re) = {0}. If however (a m , a p ) < the situation becomes more interesting. 
In the generic case (a m , a p ) = -lwe find 

Gl(re) = {a m , a m + a p }, G p (re) = {a p , a m + a p }. 

The only two exceptions of this rule for the classical series are m = r — 1, p = r for 
q ~ B r and C r . Then we have: 



G r-i(^) = {ot r -i,at r -i 
for g ~ B r and 

G^_i(k) = {a r _i,a r _i 



a r , a r _i + 2a,-}, 



a r , 2a r -i + a r }, 



G^(re) = {a r , a r -i + a r ,a r -i + 2a r }, 



G*(k) = {a r , a r -i + a r , 2a r -i + a r }, 



for g ~ C r . These last relations allow us to calculate the asymptotics of fo r all 
possible values of k for ibs = l 77 *^}. Then it is not difficult to insert them in (4.6) and 



evaluate the asymptotic behavior of all qk(t). Several examples of such calculations 
were presented in Section 4 above. 
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